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ABSTRACT 

Context. Recent investigations of the white dwarf (WD) + He star channel of Type la supernovae (SNe la) imply that this channel can 
produce SNe la with short delay times. The companion stars in this channel would survive and potentially be identifiable. 
Aims. In this Letter, we study the properties of the companion stars of this channel at the moment of SN explosion, which can be 
verified by future observations. 

Methods. According to SN la production regions of the WD + He star channel and three formation channels of WD + He star systems, 
we performed a detailed binary population synthesis study to obtain the properties of the surviving companions. 
Results. We obtained the distributions of many properties of the companion stars of this channel at the moment of SN explosion. We 
find that the surviving companion stars have a high spatial velocity (>400km/s) after SN explosion, which could be an alternative 
origin for hypervelocity stars (HVSs), especially for HVSs such as US 708. 
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Type la supernovae (SNe la) play an important role in the study 
of cosmic evolution, especially in cosmology. They have been 
applied successfully in determining cosmological parameters 
(e.g., Q. and A; Riess et al. 1998; Perlmutter et al. 1999). It is 
generally believed that SNe la are thermonuclear explosions of 
carbon-oxygen white dwarfs (CO WDs) in binaries (for the re- 
view see Nomoto et al. 1997). However, there is still no agree- 
ment on the nature of their progenitors (Hillebrandt & Niemeyer 
2000; Podsiadlowski et al. 2008; Wang et al. 2008), and no SN la 
progenitor system has been conclusively identified from before 
the explosion. 

Over the past few decades, two families of SN la progeni- 
tor models have been proposed, i.e., the double-degenerate (DD) 
and single-degenerate (SD) models. Of these two models, the SD 
model is widely accepted at present (Nomoto et al. 1984). It is 
suggested that the DD model, which involves the merger of two 
CO WDs (Iben & Tutukov 1984; Webbink 1984; Han 1998), 
likely leads to an accretion-induced collapse rather than to an 
SN la (Nomoto & Iben 1985). For the SD model, the compan- 
ion is probably a main-sequence (MS) star, a slightly evolved 
subgiant star (WD + MS channel), or a red-giant star (WD + 
RG channel) (e.g., Hachisu et al. 1996, 1999a,b; Li & van den 
Heuvel 1997; Langer et al. 2000; Han & Podsiadlowski 2004, 
2006; Chen & Li 2007, 2009; Meng et al. 2009; Lii et al. 2009; 
Wang, Li & Han 2009). An explosion following the merger of 
two WDs would leave no remnant, while the companion stat- 
in the SD model would survive and be potentially identifiable 
(Podsiadlowski 2003). There has been no conclusive proof yet 
that any individual object is the surviving companion star of an 
SN la. It will be a promising method to test SN la progenitor 
models by identifying their surviving companion stars. 



Yoon & Langer (2003) followed the evolution of a CO WD 
+ He star system with a 1 .0 M CO WD and a 1 .6 M He star 
in a 0.124d orbit. In this binary, the WD accretes He from the 
He star and grows in mass to the Chandrasekhar (Ch) mass. SNe 
la from this binary channel can neatly avoid H lines. Recently, 
Wang et al. (2009a) systematically studied the WD + He star 
channel of SNe la. In the study, they carried out binary evolution 
calculations of this channel for about 2600 close WD binaries, in 
which a CO WD accretes material from an He MS star or an He 
subgiant to increase its mass to the Ch mass. The study shows 
the parameter spaces for the progenitors of SNe la. By using a 
detailed binary population synthesis (BPS) approach, Wang et al. 
(2009b) find that the Galactic SN la birthrate from this channel 
is ~0.3 x 10" 3 yr _I and that this channel can produce SNe la with 
short delay times (~45-140Myr) from the star formation to SN 
explosion. The companion star in this channel would survive and 
show distinguishing properties. 

In recent years hypervelocity stars (HVSs) have been ob- 
served in the halo of the Galaxy. HVSs are stars with a velocity 
so great that they are able to escape the gravitational pull of the 
Galaxy. However, the formation of HVSs is still unclear (for a 
recent review see Tutukov & Fedorova 2009). It has been sug- 
gested that such HVSs can be formed by the tidal disruption of 
a binary through interaction with the super-massive black hole 
(SMBH) at the Galactic center (GC) (Hills 1988; Yu & Tremaine 
2003). The first three HVSs have only recently been discov- 
ered serendipitously (e.g., Brown et al. 2005; Hirsch et al. 2005; 
Edelmann et al. 2005). Up to now, about 17 HVSs have been 
discovered in the Galaxy (Brown et al. 2009; Tillich et al. 2009), 
most of which are B-type stars, probably with masses ranging 
from 3 to 5 M (Brown et al. 2005, 2009; Edelmann et al. 2005). 
One HVS, HE 0437-5439, is known to be an apparently normal 
early B-type star. Edelmann et al. (2005) suggests that the star 
could have originated in the Large Magellanic Cloud, because it 
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is much closer to this galaxy (18 kpc) than to the GC (see also 
Przybilla et al. 2008). 

At present, only one HVS, US 708, is a subdwarf O (sdO) 
star, and Hirsch et al. (2005) speculatS that US 708 is formed by 
the merger of two He WDs in a close binary induced by the in- 
teraction with the SMBH in the GC and then escaped. Recently, 
Perets (2009) has suggested that US 708 may have been ejected 
as a binary from a triple disruption by the SMBH, which later 
on evolved and merged to form an sdO star. Because of the ex- 
istence of the short orbital periods (~1 h) for the WD + He star 
systems, Justham et al. (2009) argues that the WD + He star 
channel of SNe la may provide a natural explanation for stars 
like US 708. 

Han (2008a) obtained the distributions of many properties of 
the surviving companions from the WD + MS channel of SNe 
la. The properties can be verified by future observations. The 
purpose of this Letter is to investigate the properties of the sur- 
viving companions of the WD + He star channel and to explore 
whether HVSs such as US 708 could have been released from 
the binaries that produced SNe la. In Section 2, we describe the 
BPS approach and the simulation results for the properties of the 
surviving companions. Finally, a discussion is given in Section 
3. 



2. The distributions of properties of the companion 
stars 

In the WD + He star channel, the progenitor of an SN la is a 
close WD binary system, which has most likely emerged from 
the CE evolution (Paczyhski 1976) of a giant binary system. The 
CE ejection is still an open problem. Here, we use the standard 
energy equations (Webbink 1984) to calculate the output of the 
CE phase. The CE is ejected if 



2fl f 



2a i 



^don 



(1) 



where A is a structure parameter that depends on the evolution- 
ary stage of the donor, M^ on is the mass of the donor, M acc the 
mass of the accretor, a the orbital separation, M env the mass of 
the donor's envelope, Rd on the radius of the donor, and the in- 
dices i and f denote the initial and final values, respectively. The 
right side of the equation represents the binding energy of the 
CE, the left side shows the difference between the final and ini- 
tial orbital energy, and ff ce is the CE ejection efficiency, i.e., the 
fraction of the released orbital energy used to eject the CE. For 
this prescription of the CE ejection, there are two highly uncer- 
tain parameters (i.e., A and ff ce ). As in previous studies, we com- 
bine a ce and A into one free parameter a ce A, and set it to be 0.5 
(e.g., Wang et al. 2009b). 

To obtain the distributions of properties of the surviving 
companions, we performed a Monte Carlo simulation in the BPS 
study. In the simulation, by using the Hurley's rapid binary evo- 
lution code (Hurley et al. 2000, 2002), we followed the evolution 
of 4 x 10 7 sample binaries from the star formation to the forma- 
tion of the WD + He star systems according to three evolutionary 
channels (i.e., the He star channel, the EAGB channel, and the 
TPAGB channel; for details see Wang et al. 2009b). If a binary 
system evolves to a WD + He star system, and if the system, 
at the beginning of the Roche lobe overflow (RLOF) phase, is 
located in the SN la production regions in the plane of (logP 1 , 
Mj,) for its M\^ D , where P\ M\ and M^, D are, respectively, the 
orbital period, the secondary's mass, and the WD's mass of the 
WD + He star system at the beginning of the RLOF (see Fig. 
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Fig. 1. The distribution of properties of companion stars in the 
plane of (V rb, M| N ) at the current epoch, where V or b is the or- 
bital velocity and M? N the mass at the moment of SN explosion. 



8 of Wang et al. 2009a), we assume that an SN la is produced, 
and the properties of the WD binary at the moment of SN explo- 
sion are obtained by interpolation in the three-dimensional grid 
(M^, D , M\, logP') of the -2600 close WD binaries calculated in 
Wang et al. (2009a). 

In the BPS study, the primordial binary samples are gener- 
ated in the Monte Carlo way. We adopted the following input for 
the simulation (e.g., Han et al. 1995a, 2002, 2003, 2007; Han 
2008b; Wang et al. 2009b). (1) The star-formation rate (SFR) 
is taken to be constant over the past 15Gyr. For the constant 
SFR, we assume that a binary with its primary more massive 
than 0.8 M is formed annually (e.g., Han et al. 1995b). (2) The 
initial mass function (IMF) of Miller & Scalo (1979) is adopted. 
(3) The mass-ratio distribution is taken to be constant. (4) The 
distribution of separations is taken to be constant in log a for 
wide binaries, where a is the orbital separation. (5) The orbits 
are assumed to be circular. 

The simulation gives current-epoch distributions of many 
properties of companions at the moment of SN explosion, e.g., 
the masses, the orbital periods, the orbital separations, the orbital 
velocities, the effective temperatures, the luminosities, the sur- 
face gravities, the surface abundances, the mass-transfer rates, 
the mass-loss rates of the optically thick stellar winds, etc. The 
simulation also shows the initial parameters of the primordial 
binaries and the WD binaries that lead to SNe la. Figures 1-5 
are selected distributions that may be helpful for identifying the 
surviving companion stars. 



3. Discussion 

Figures 1 and 2 show the distributions of the masses, the orbital 
velocities, the effective temperatures, and the surface gravities 
of companion stars at the moment of SN explosion. In Figure 
1, the companion star has an orbital velocity of ~300-500km/s 
for a corresponding mass of ~0.6-1.7M o at the moment of SN 
explosion. In the three formation channels of the WD + He star 
systems, SNe la mainly come from the He star channel. For this 
formation channel, the recorded properties at each step show that 
a primordial binary system with a primary mass M\\ ~ 5.0 - 
8.0 M e , a secondary mass M-^ ~ 2.0 - 6.5 M e , and an orbital 
period P; ~ 10 - 40 d would evolve to a close WD + He star 
system with a WD mass M' ~ 0.87 - 1.2 M G , a He star mass 



Ml 



1.0 - 2.6 M Q , and an orbital period P 1 ~ 0.04 - 0.2 d. 



B. Wang & Z. Han: Companion stars of SNe la and HVSs 



3 




5.3 5.2 5.1 5 4.9 4.8 4.7 4.6 

iog(r eft /K) 

Fig. 2. Similar to Fig. 1, but in the plane of (log T e s, log g), where 
T e g is the effective temperature of companion stars at the mo- 
ment of SN explosion and logg the surface gravity. The com- 
panion stars are out of thermal equilibrium, e.g., a 1.243 M 
He star with 0.274fl o and (logoff, logg)=(4.70, 5.65) will be 
a 0.267 R He star and with (log r eff , logg)=(4.75, 5.67) after 
the He star back to the thermal equilibrium. 
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Fig. 3. The distribution of spatial velocity for surviving compan- 
ion stars of SNe la. The dotted line denotes the results of ejecta 
velocity HOOOkm/s, while the solid line shows ejecta velocity 
13500km/s. 
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Fig. 4. Similar to Fig. 1, but in the plane of (logf SN , M| N ), 
where f SN is the orbital period at the moment of SN explosion. 
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Fig. 5. Similar to Fig. 1, but in the plane of (V rot , M| N ), where 
V rot is the equatorial rotational velocity of companion stars at the 
moment of SN explosion. 

Finally, the WD + He star system results in an SN la explosion 
and survives a companion star. 

However, Figures 1 and 2 are for the moment of SN explo- 
sion, which could be modified by the explosion. The SN ejecta 
will interact with its companion. The companion stars will be 
stripped of some mass and receive a kick velocity that is per- 
pendicular to the orbital velocity. Adopting a similar method 
to Meng, Chen & Han (2007), we estimated the stripped mass 
for the companion stars of the WD + He star channel and find 
that the stripped mass is very low, e.g., a 1.243 M companion 
star with 0.274 R G only loses a mass of 0.015 M in our simula- 
tion. We also roughly obtain the kick velocity based on the mo- 
mentum conservation equation^ The kick velocity mainly de- 
pends on the ratio of separation to the radius of companions 
at the moment of SN explosion, A/R^ N , and the leading head 
velocity of SN ejecta, which is assumed to be in the range of 
11000 to 13500km/s. The velocity of HOOOkm/s is from the 
SN ejecta kinetic energy 1.0 x 10 51 erg corresponding to the 
lower limit of normal SN la kinetic energy, while the velocity 
of 13500km/s is from the SN ejecta kinetic energy 1.5 x 10 51 erg 
corresponding to the upper limit of the kinetic energy (Gamezo 
et al. 2003). We can obtain the spatial velocity by the formula 

V 2 N = V^ick + ^orb' where ^kick and V orb are the kick velocity 
and the orbital velocity of the companion star at the moment of 
SN explosion, respectively. In Figure 3, we show the distribu- 
tion of the spatial velocity for the surviving companion stars of 
the WD + He star channel. We see that the surviving companion 
stars have high spatial velocities (>400km/s) that almost entirely 
exceed the gravitational pull of the Galaxy nearby the sun. Thus, 
the surviving companion stars from the WD + He star channel 
could be an alternative origin for HVSs. 

US 708 is an extremely He-rich sdO star in the Galaxy halo, 
with a heliocentric radial velocity of +708 + 15 km/s (Hirsch et 
al. 2005). We note that the local velocity relative to the Galatic 
center may lead to a higher observation velocity for the sur- 
viving companion stars, but this may also lead to a lower ob- 
servation velocity. Considering the local velocity near the sun 
(~220km/s), we find that ~30% of the surviving companion 
stars may be observed to have velocity V > 700 km/s for a given 



1 However, in reality, the collision between the SN ejecta and its 
companion cannot be elastic. If we adopt the inelastic collision in our 
calculations, the kick velocity will only decrease by ~3%-5%. Thus, 
the inelastic collision is no significant influence on the kick velocity. 
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SN ejecta velocity 13500km/s. In addition, the asymmetric ex- 
plosion of SNe la may also enhance the velocity of the surviv- 
ing companions. Thus, a surviving companion star in the WD + 
He star channel may have a high velocity like US 708 (see also 
Justham et al. 2008). 

The companion stars are out of thermal equilibrium at the 
moment of SN explosion. For He stars, the equilibrium radii 
are lower than at the moment of SN explosion. Thus, the sur- 
face gravity at equilibrium should be greater than the one in 
Figure 2; e.g., a 1.243 M e He star with 0.274 R Q and (logr eff , 
logg)=(4.70, 5.65) will be a 0.267 R Q He star and with (log T eB , 
logg)=(4.75, 5.67) after the He star is back to the thermal equi- 
librium. It is well known that a shock will develop after impact 
by the ejecta. A large part of the material in the companion's 
envelope is heated by the shock, and some of the material is va- 
porized from the surface of the companion stars. Following a 
similar method to Chen & Li (2007), we also estimated the va- 
porized mass, and find that the mass loss from the companion 
stars is not significant (<5%). Such a surviving companion star 
may be significantly overluminous or underluminous depending 
on the amount of heating (e.g., Podsiadlowski 2003). Figure 2 
could be a starting point for further studies of this kind. 

Figure 4 shows the distributions of orbital periods and sec- 
ondary masses of the WD + He star systems at the moment of 
SN explosion. The orbital periods and secondary masses of the 
WD + He star systems at this moment are basic input parame- 
ters when one simulates the interaction between SN ejecta and 
its companion. It is suggested that, for hot stars with radiative 
envelopes (such as He stars), tidal forces may be inefficient for 
synchronization (e.g., Zahn 1977). However, the recent study by 
Charpinet et al. (2008) supports efficient tidal synchronization 
for hot subdwarf stars. The work by Toledano et al. (2007) also 
indicates, on the other hand, that even stars with radiative en- 
velopes may have efficient tidal interaction on a time scale com- 
parable to convective envelopes. Thus, we make an assumption 
that the companion stars co-rotate with their orbits. In Figure 5, 
we show the distributions of equatorial rotational velocities of 
the companions stars. We see that the surviving companion stars 
are fast rotators, so their spectral lines should be broadened no- 
ticeably. 

The simulation in this Letter was made with a ce A = 0.5. 
If we adopt a higher value for a ce A, e.g., 1.5, the birthrate of 
SNe la would be a little bit higher and the delay time from the 
star formation to SN explosion longer. This is because binaries 
emerging from CE ejections tend to have longer orbital periods 
for a large a ce A and are more likely to be located in the SN la 
production region (Fig. 8 of Wang et al. 2009a). Due to the lack 
of WD binaries with short orbital periods (logP 1 < -1.2) for a 
large a ce A, the companions with orbital velocity V or b > 430 km/s 
would be noticeably absent in Figure 1 . 

The distributions are results of the current epoch for a con- 
stant SFR. For a single starburst, most of the SN explosions oc- 
cur between ~45 Myr and ~ 140 Myr after the starburst; i.e., SNe 
la from the WD + He star channel will be absent in old galaxies. 
The Galactic SN la birthrate from this channel is ~0.3xl0~ 3 yr 1 
(Wang et al. 2009b). By multiplying the birthrate with a typical 
MS lifetime of He stars, ~10 7 yr, we estimated the current num- 
ber of this type of HVSs in the Galaxy to be ~10 3 . In future 
investigations, we will employ the Large sky Area Multi-Object 
fiber Spectral Telescope (LAMOST) to search the HVSs origi- 
nating from the surviving companion stars of SNe la. 
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